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Abstract. Wc construct a regularized cumulant (semi-invariant) representation of a solution 
of the initial value problem for the BBGKY hierarchy for a one-dimensional infinite system 
of hard spheres interacting via a short-range potential. An existence theorem is proved for 
the initial data from the space of sequences of bounded functions. 
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1 Introduction 

Last ten years new functional-analytical methods of investigation of the initial value problem to 
the BBGKY hierarchy for infinite systems have been developed in the books 0121121 13 • As is 
well known when constructing a solution of the BBGKY hierarchy of such systems, considerable 
analytical difficulties arise [31 IHl 13 IH1 • Until recently, existence theorems have been proved for 
one-dimensional systems in the case of a certain class of short-range interaction potentials and 
three-dimensional case only for hard spheres systems 0. 

While constructing the solution of the initial value problem for the BBGKY hierarchy of 
the classical systems of particles with the initial data from the space of sequences of bounded 
functions, one is faced with certain difficulties related to divergence of integrals with respect to 
configuration variables in each term of an expansion of the solution [5] (see also 12 ) . The 
same problem arises also in the case of the cumulant representation of the solution stated in [2] . 

In this paper, we propose a regularization method for the solution of the BBGKY hierarchy 
in the cumulant representation. Due to this method, the structure of the solution expansions 
guarantees the mutual compensation of the divergent integrals in every term of the series. We 
establish convergence conditions for the series of the solution and prove an existence theorem 
of a local in time weak solution of the BBGKY hierarchy for the initial data from the space 
of sequences of functions which are bounded with respect to the configuration variables and 
exponentially decreasing with respect to the momentum variables. 

2 Initial value problem for BBGKY hierarchy 

Let us consider a one-dimensional system of identical particles (intervals with length a and unit 
mass m = 1) interacting as hard spheres via a short range pair potential Every particle i is 
characterized by phase coordinates {qi,Pi) = € M x M, i > 1. For the configurations qi € of 
such a system (qi is the position of the center of the iih particle), the following inequalities must 
be satisfied: \qi - qj\ > a, i ^ j > 1. The set Wn = {{qi, ■■■ ,qn} \ 3(i, j), i ^ j e {1, . . . ,n} : 
IQi ~ Qj\ < defines the set of forbidden configurations in the phase space of a system of n 
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particles. The phase trajectories of such hard sphere system are determined almost everywhere 
in the phase space {xi, . . . ,Xn} € M" x (R" \ Wn), namely, outside a certain set A^^ of the 
Lebesgue measure zero p. The initial data {xi, . . . , Xn} G M" x (M" \ Wn) belong to the set 
if (a) there is more than one pair collision at the same moment of time t € (—00, +00) or 
(b) infinitely many collisions occur within a finite time interval. 

We assume that the interaction between the hard spheres is given by a potential $ with 
a finite range R such that the following conditions are satisfied: 



(a) $ G C^{[a, R]), < a < R < 00, 

+00, \q\ € [0, cr), 

^ 0, \q\ € (i?,cx)), 

(c) $'((j + 0) = 0. (1) 



(b) mq\) 



We note that conditions ^ imply the estimate 



i<j=i 



<bn, sup I |([f]), (2) 

qe[a,R] 



where [^] is the integer part of the number ^. 

The evolution of states of the system under consideration is described by the initial value 
problem for the BBGKY hierarchy ^ |2j 

d 
ot 



, Xl , . . . , Xg-^-l j 

s °? 

+ / dPP ■ (^Fg^i{t,xi, ...,Xg,qi + a, Pi - P) 



i=l Q 



- Fg+i{t,xi, . . . ,qi,pi - P, . . . ,Xg,qi - a,pi) 
+ Fg+i{t,xi, ...,Xg,qi- a, Pi + P) 

- Fg+i{t,xi, . . .,qi,Pi + P, . . . ,Xg,qi + cr,pij^, (3) 
Fg{t,Xl,...,Xs)\^^^^ = Fs{0,xi,...,Xs), s > I, (4) 

where {•, •} is the Poisson bracket, Hg is the Hamiltonian of the s particle system, and F{Q) = 

),...) is a sequence of initial s-particle distribution functions 
Fg{0,xi, . . . , Xg) defined on the phase space W x {W \ Wg). 

Consider the initial value problem for the BBGKY hierarchy ^ with the initial data F{0) 
from the space L|°^ of sequences / = (l, . . . , fn{xi, ■ ■ ■ , . . . ) of bounded functions 

fn{xi, . . . , Xn), /o = 1, > 0, that are defined on the phase space M"" x (M" \ Wn)-, are invariant 
under permutations of the arguments Xi, i = 1, . . . ,n, and are equal to zero on the set Wn- The 
norm in this space is defined by the formula 



SUpC " sup \fn{xi,...,Xn)\exp\pS2^\^ (5) 
n>0 xi,...,x„ ' '-^ ' 
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where ^, (3 are positive integers. Note that the sequences of n particle equilibrium distribution 
functions of infinite systems belong to the space [Tl llUj. 

Let Y = {xi, . . .,Xs}, X = {Y,Xs+i, ■ ■ ■ ,Xs+n}, namely, X \ Y = {xg+i, ■ ■ ■ ,Xs+n}, and let 
the symbol \X\ = \Y\ + \X \ Y\ = s + n denote the number of the elements of the set X. 
By the symbol Xy we denote the set X with the subset Y treated as a single element similar 
to Xs+1t • • • 1 Xs+n- Let Lq be the subspace of finite sequences of continuously differentiable 
functions with compact supports of space of sequences of integrable functions. For -F(O) 
from Lq, and hence for all -F(O) G Lq Pi L'^^ it was proved in [SI IH] that the solution F{t) = 
(l, Fi{t, xi), . . . , Fs{t, xi, . . . , Xs), . . . ) of the initial value problem (jSJ, @ is determined by the 
series expansion 

oo 

F\Y\{t,Y) = Y,^^ J d{X\Y) (-1)'""'(|P|-1)! 

R"x(R"\H^„) P:XY=[jXi 

^ U S\x,\{-t,Xi)F\x\{0,X), \X\Y\>0, (6) 

XiCF 

where Yl is the sum over all possible partitions P of the set Xy into |P| mutually disjoint 
p 

nonempty subsets Xi C Xy, X^ D Xi = 0, k ^ I, such that the entire set Y is contained in one 
of the subsets Xi. 

On the set of sequences / G LqPiL^^, the evolution operator S\Xi\{—t, Xi) from expansion © 
is given by the formula 



iS{-t)f)^x,iiXi) = {Six,i-t)fix,\XXi) ^ Six,\{-t,Xi)f^x,\{Xi) 

'f\x,i{M-t,Xi),...,X^x,\{-t,Xi)), ifxG(Rl^'lx (Ml^'l\VF|x,|))\XO,^|, 



[O, if xGRI^'Ix 

where Xi{—t,Xi), i = 1, . . . , l-'^^zl, is the solution of the initial value problem for the Hamilton 
equations of the system of \Xi\ particles with initial data Xi(0,X;) = xi (5|X;|(0) = / is the 
identity operator). Under conditions on the potential the evolution operator ^ exists 
for t G (— oo, +oo); its properties are described in p. 

In the nth term of expansion @, the form of the integrands is constructed by using the 
cumulant of order 1 + n for the evolution operators ((TJ: 

\{S\x,\{-t,Xi)^^^^\x\Y\{t,XY), |x\y|>o. (8) 

P:Xy=UX, XiCP 

I 

Here, the notation from formula © is used. Note that the order of the cumulant ^i+\x\Y\{t) is 
determined by the number of elements of the set Xy (in this case, by 1 + |X \ y| elements). 

The simplest examples (1 + |X \ y| = 1, 2, 3) of the evolution operator 2li+„(t) Q have the 
form linn HI] 

%i{t,Y) = Ss{-t,Y), (8a) 
^2{t,Y,x,+i) = Ss+i{-t,Y,x,+i) - Ss{-t,Y) Si{-t,xs+i), (8b) 

"^zit^Y, Xs+l,Xs+2) = Ss+l{-t,Y,Xs+l,Xs+2) " Ss+l{-t,Y,Xs+l) Si{-t,Xs+2) 
- Ss+l{-t,Y,Xs+2) Sl{-t,Xs+l) - Ssi-t,Y) S2{-t,Xs+l,Xs+2) + 

+ 2\Ss{-t, Y) Si{-t, Xs+i) Si{-t, Xs+2). (8c) 

Thus, the cumulant representation of the solution @ of the initial value problem for the BBGKY 
hierarchy ©, (0]) is determined by the cumulants 2ti_|_|x\y|(i) (© for the evolution operators 
S\x,\i-t,Xi) 0. 
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3 Regularization of solution 

For -F(O) E every term n = \X \Y\ of expansion © contains divergent integrals with 

respect to the configuration variables. Let us show that the above-stated cumulant nature 
of the solution expansions © for the initial value problem of the BBGKY hierarchy Q, @ 
guarantees the compensation of the divergent integrals, i.e., the cumulants are determined terms 
of expansion ® as the sum of summands with divergent integrals that compensate one another. 
In order to prove this fact, let us rearrange the terms of expansion so that they are represented 
by the simplest mutually compensating groups of summands. Such procedure will be called 
a regularization of the solution @. In this case, the regularization will be based on expressing 
cumulants of higher order in terms of the first and second order cumulants. For fixed initial 
data the second-order cumulants will be determined by expressions which compensate each other 
over a certain bounded domain. 

The next lemma shows that, in the general case, the cumulant of the (1 -|-n)th order 
n > 1, is expressed via the first and second order cumulants. 

Lemma 1. The equality 

^i+lx\Y\{t,XY)= Yl Mt,Y,Z) 

ZCX\Y 

X E 

P:X\{YUZ)=\J X, XiCP 



^ (-1)'"(|P|)! n 2ti(t,^/), I^\>^I>1, (9) 



is true, where ^ is the sum over all the nonempty subsets Z of the set X \Y, Z C X \Y, 
z 

and the group of \Z\ particles evolves as a single element, and ^ is the sum over all possible 

p 

partitions P of the set X\{Y L) Z) into |P| mutually disjoint nonempty subsets Xi C X\{Y U Z), 
Xk Xi = 0, k ^ I, such that every cluster of \Xi\ particles evolves as a single element. 

The proof of Lemma ^ is based on the verification that © is equal to expression by 
taking into account the representation of a second order cumulant in terms of the first order 
ones. 

By using Lemma ^ represent integrands of every summand from expansion in terms of 
the first and second order cumulants. As a result, expansion © for the solution of the initial 
value problem ©, @ takes the form 

oo 

=2ti(t,y)F|y|(o,y) + ^- / d{x\Y) (lo) 



'""-^ E"x(R"\W„) 



X 



Mt,Y,z) Y {-if^\p\\ llMt,Xi)F\xi{0,X), |X\y|>L 



ZCX\Y P:X\(YUZ)=\J Xi XiCP 

Zj^0 I 

For the initial data F(0) G n L|°^, the equality 



/ 



diX\Y) Y Mt,Y,Z) Y (-1)'"|P|! n 2ti(i,^;)i^|x|(0,X) 

w ) ZCX\Y P:X\(YUZ)=[JXi XiCP 

J diX\Y) Y (-l)l''\(''^^%2(t,>^,^)i^|x|(0,X), \X\Y\>1, (11) 



^ Z+0 
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is true. Here, we have used the Liouville theorem and taken into account the relation 

m 

s(m, k) = (-1)™, m > 1, (12) 

k=l 

where s{m, k) is the Stirling number of the second kind defined as the number of all distinct 
partitions of a set containing m elements into k subsets. 

Thus, by virtue of equality the expansion (jlUj) for the solution of the initial value 

problem for the BBGKY hierarchy (jSJ, @ takes the form 

F|y|(t,y) =2ti(t,y)F|y|(o,y) (is) 

oo 

+ J d{X\Y) (-l)l''\("'^^%2(i,>^,^)i^x|(0,X), \X\Y\>1, 

where the notation from Q is used. 

Representation (|13() will be called a regularized cumulant representation of the solution of 
the initial value problem for the BBGKY hierarchy (|3|). 



4 Existence theorem 

Taking into account the invariance of the Gibbs distributions (the Maxwellian distribution can 
be extended to become a Gibbs distribution) with respect to the action of evolution opera- 
tors S (— t) and using the condition Q and the relation Q we estimate the integrands in 
expansion (jlSf) . 

Lemma 2. If F{0) G then the inequality 



(-l)l^\(''^^%2(t,l^,^)i^x|(0,X) 

ZCX\Y 
Z^0 




< 2||F(0)|| (^e2^^)^(ee2^^)"exp^ -/3y ^ I (14) 



holds, where the notation from formulae ((2)), @ and Q is used. 

As a consequence of Lemma 121 the following existence theorem is true. 

Theorem. // F[Q) G is a sequence of nonnegative functions, then for ^ < ^~^~~\f^ 

and t G [0,to), where to = ^ ' \/^ ~ ^0 ' ^ max(2i?, 1), C2 = max (2(46 + 1), 

P = P' + P" , b = sup 1*^(9)1 ([■§]) [^] is the integer part of the number ^, there exists 

qe[a,R] 

a unique weak solution of the initial value problem for the BBGKY hierarchy ©, Q, namely, 
the sequence F{t) € L|°^ of nonnegative functions Fs{t) determined by expansion (jlSp . 

Proof. Let particles interact via a short-range pair potential that satisfies conditions We 
assume that, at the initial instant, the configuration coordinates gj, i = 1, 2, . . . , s, of the particles 
constituting the cluster Y take values in a compact set of those |y| intervals k with length |/y| 
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that Qi S /j. Then if during the time interval [0, t) none of the particles of any cluster Z G X\Y 
interacts with the particles of cluster Y, then the operator equality holds 

SiYuzi i-t, Y, Z) = S\y\ i-t, Y)S\zi 
and, as a result, we have 

2i2(t,y,z)i^x|(o,x) = 0. 

Therefore, in this case, the integrands in the nth term of expansion (jl^j) are equal to zero. Since 
Z C X\Y, it follows that in the expansion (|13|) the domain M" \ Wn of integration with respect 
to the configuration variables is determined by n bounded intervals where the particles of the 
cluster Y during the time interval [0, t) interact with particles of the cluster Z that contains the 
maximal number of particles, namely, with the particles of the cluster X\Y. Thus, the domain 
M" \ Wn of integration estimates the following finite value 

(s+n \ " 
C + Cot + {Ci + C2t)n + t ^ , (15) 

i=s+l J 

s 

In view of estimates (fTl)) and (fT3|l , we obtain 



where C =\Iy\ + 2sR, Cq = 2s(46 + 1) + E pI Ci = 2R, and C2 = 2(46 + 1). 

i=l 



i=l J n=0 

X / dps+i ■ ■ ■ dps+n expl-P^^\iC + Cot + (Ci + C2t)n + t ^ 

-'^'^ I i=s+l J V i=s+l 



Pi 



Taking in account the relation 

(s+n 
{C + Cot) + (Ci + C2t)n + t E 
i=s+\ 

k=0 ' r=0 ' ^ \i=s+l 

and the inequality 

s+n \ I 2 1 /9\' 

we compute the integrals with respect to the momentum variables in the right-hand side of ()16() : 

/C s+n 2 ^ / 2tt\^ 

dps+i ■ ■ ■ dps+nexp I - /5" E y I = (^J ' /9 = /5' + f3". (18) 

i=S + l 

Then estimate H16|) takes the form 



s+n \nkr ^ ^ ^ . s n-k- 



I (t, Y) I < 2 ||F(0) II {Ce'f^'r exp | -/3 E f 
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n=0 ^'^ ' k=0 r=0 ^ ' 



Let us put Ci = max(Ci, 1) and C2 = max ^6*2, Jrj and continue estimate (^fl). For arbitrary 
t > 0, the inequalities Ci + 6*2^ > 1 and {Ci + C2t)^ > 1 are true, and, therefore. 



By using the inequalities ^ e'^ and X] < e'-'-^"'"'-^''*^ and taking estimate (|16l|) into 

r=o ■ fc=0 

account, we obtain 



< 2||F(0)|| (^e2^'')«exp 



El. 

2 



(C+Cot) ^ j 2ee2/^^+iy^ ] (Ci + C2i)"- (19) 



X e 

n=0 



Thus, if ^ < y then series ^ converges for < t < to = ^ (^^^ V Cf " • 

We have thus shown that, under the conditions of the theorem, series 1)16^ converges. 

Finally, using the well-known theorems of functional analysis |13j and arguing similarly 
to we show that the sequence F{t) is the unique weak solution of the initial value problem 

for the BBGKY hierarchy ©, ©. ■ 



5 Conclusion 

The theory of the BBGKY hierarchy is developing now since the area of its application grows jl4l 
I15j . At present while solving the initial value problem of the BBGKY hierarchy, many math- 
ematical problems arise jSl El El ^1- In this paper one of such problems is considered for 
infinite particle system employing the method of the interaction region . Taking into account 
the cumulant representation [211^1 we construct a new regularized representation of the solution 
of the BBGKY hierarchy for a one-dimensional infinite system of hard spheres interacting via 
a short-range potential. For the initial data from the space of sequences of functions which 
are bounded in configuration variables and exponentially decreasing in momentum variables, 
existence of the cumulant representation allows to properly regularize such expansion for the so- 
lution, determining by second order cumulant as the sum of summands with divergent integrals 
that compensate one another. A multidimensional case and a case of a more general interaction 
potential will be investigated in further contribution. 
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